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Abstract 



43 

. Thermodynamic stable interaction pair potentials which are not of the 

form "positive function + real continuous function of positive type" are 
presented in dimension one. Construction of such a potential in dimension 
two is sketched. These constructions use only elementary calculations. The 
mathematical background is discussed separately. 

> 

^ ■ PACS numbers: 05.20.-y, 02.20.-a, 02.40.Ft 

o, 

f*"^ , Keywords: thermodynamic stability, convex cone 

o 

> 

• rH . 

^ ■ 1 Introduction 

. P.. 

In Ruelle's book |R69j on statistical mechanics, in section 3.2 concerning one 
species of classical particles in W, you can read: 

1 PROPOSITION. // the pair potential $ can be written in the form 

$ = $ 1 + $ 2 (1) 

where $i is positive, and $ 2 is a real continuous function of positive type, then 
$ is stable. 
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"Positive" is meant here and throughout this paper as nowhere negative, 
"stable" means 

3E eR such that VN, V{ Xl ...x N } C W : U(xf-x N ) ^ N ■ E , (2) 

where U(x± ■ • • xjy) = ®{ x j ~ x i)- (3) 

i¥=3 

This proposition is accompanied by the 

2 FOOTNOTE. It seems to be an open problem to construct a stable potential 
which is not of the form (QJj- 

We solve this problem in dimension 1, considering particles either in Z or in 
R, giving a detailed proof. In dimension 2 the problem can also be solved, but 
we give only a sketch of the ideas. 

To make it simple, we consider only pair potentials which are bounded con- 
tinuous functions and state the stability property as 

3 DEFINITION. A bounded continuous real valued function V on W is stable, 
if 

E(p) := J J p(x)V(x - y)p{y)d v xd v y > (4) 

for every positive finite measure p{x)d v x onW ' . A bounded real valued function 
V on Z u is stable, if 

E(p) :=^^p(rh)V{m-n)p(n) >0 (5) 

m n 

for every positive bounded function p(m) on 

The stability property used in Ruelle's Theorem is an immediate consequence. 
With p = J2 i=1 S(xi — Xj) put into equation (jlj) one gets 

U( Xl ■■■x N ) = E(p) - N ■ 7(0) > —N ■ 7(0). 

The main result of our considerations is stated as 

4 THEOREM. Each of the following functions is a stable pair potential, but 
not a sum of a positive and a real valued positive definite function. 

1. The function V : Z — > IR, defined as 

7(0) = 7(2) = 7(-2) = 1, 7(1) = 7(-l) = -1, (6) 
7(n) = Vn with \n\ > 3, 

2. The function W : IR — > R, defined as 

/oo 
V(n)f(n-x + y)f(y)dy, (7) 
-oo 



with f a positive continuous function (—2,0)^^ an d V as defined in (0|). 



1 Construction in higher dimensions is still an open problem. 
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2 Properties of the interaction potentials 

Proof. Of part (QJ) o/[|] Theorem. 

Denote the distribution of particles on the chain by the "density" p, a function 
Z — > Z + . The interaction energy U becomes smaller, when the system is cut into 
non-interacting pieces: If p(n) > p(n + 1) divide the chain, cutting between n + 1 
and n + 2. Moving the pieces apart, one looses the energy 

2[p(n) - p{n + l))p{n + 2) + 2p(n + l)p(n + 3) > 0. 

The symmetric procedure of cutting between n — 2 and n — 1 lowers the energy 
if p(n — 1) < p{n). 

Now there remains a set of pieces of no more than three lattice points, with 
densities like 

< p(n - 1) < p{n) > p{n + 1) > 0. 
Including the "self-energies" iV • 1^(0) one gets for each piece, centered around n, 

E = p(n-l) 2 + p(n) 2 + p(n + l) 2 + 2[p(n-l)p(n + l)-p(n-l)p(n)-p(n)p(n + l)} 

= [p{n - 1) - p(n) + p(n + l)] 2 > 0. 

Proving the stability of V. 

If V were the sum of a positive and a positive definite function, it would give 

5>(nHn)>0, (8) 
for each p being both positive and positive definite. Now consider 

/i(5z/) = l, p(5v±l) = y^y±, M (5z/±2) = 0, (9) 

which is obviously positive. Positive definiteness is seen by using Bochner's the- 
orem [RN55] and calculating the Fourier- Transform, with a G (— 7T, 



n 

2tt 

T 



n)e~ ma 



V55(a) + [5{a - —), +6(a + — ) J 



>0. (10) 



But it does not give a positive value in (jSJ): 

^V(n)/i(n) = 2 - v 7 ^ < 0. 

n 

□ 
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The appearance of the numbers 5 and a/5 may seem mysterious. Demystifying 
is the next section, where we present the "origin" of these V and p. 

In this section we develop further use of these functions in R and in R 2 . 

Proof. Of part (TJj] of\4\ Theorem. 

For iV particles at xi . . . x^ consider the measure 

!>(■'■) ^(.r -,-,). (11) 

j 

Adding the self-energies N • W(0), we study 

E — J J p(x)W(x — y)p(y) dxdy 

= J2v(n) f p f (x + n)pf{x)dx, (12) 

n J 

with pf(x) := f f(x — y)p(y)dy. Splitting the integral in (fT2|) into pieces of 
intervals with unit length and defining p f tX (m) = pf(x + m) gives 

E = / dx V{n)pf(x + m + n)pf(x + m) 

meZ n 

= dx piM v (p - m )pfA m ) > °> 

■'0 ^ mgZ 



by part ([TJ) of the theorem. So the potential W is stable. 
Now consider the distribution 

P<d{x) = p,(m)5(x — m), (13) 

m 

using the sequence p defined in (jUJ). This distribution is positive and positive 
definite, as can be seen at its Fourier transform, which is (up to a factor) the 
same as in (ITU]) , now with pD(a + 2n) = pd(&) periodically extended to all a G R. 
This pd is used to show that the potential is not a sum of positive and positive 
definite functions: 

J W(x)po(x)dx 

= yV(n)y%(m) / dy f dx 5{x - m)f(n - x + y)f(y) 

™ J-l Jm. 



Y,V{n)p{n) ■ ! f 2 (y)dy < 0. (14) 

n J 

In the last step the final support of / is essential. □ 
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Construction of a stable pair potential in M 2 being a function of the particle 
distances only may be done in the following way: 

• Use W(x) defined in ([7j), now with an / supported on (—7, 7), convolute it 
twice with the distribution 



h(x) = ^V e|n| 5(x - 5n) 



W 1 {x)= / h{x-y)W{y- z)h{z)dydz. 




• Take the mean value (times 2n) of all rotated versions: W r (x) = -Wi(\x\). 

• Smoothen out W r with a positive continuous function g(r) with support on 

W 2 (x)= f J g{\2-m Wx $J^\ {\A)*V**- 

That the stability is not destroyed by the double convolution with h follows 
from a consideration as it is used in the equation (Tl2|) . Written in a formal way: 

(p\ Wi \p) = (p\h*W*h^ \p) = (p*h\W\p*h). 

Considering only smooth densities p(x) one may take Wi(xi)5(x2) as a stable 
distribution in M 2 : 

(p\ Wl ■ 5 \p)dim=2 = J (Py \ Wi \p y )dim=ldy > 0. 

Now rotating the axes and taking the mean value does not destroy the stability. 
Once more a double convolution is done, now with g in order to get W2 as a 
bounded continuous potential acting in M 2 . 

(P\ w 2 |P> = (p\ 9 * W * g- |p) = (p * g\ W \p * g) > 0. 

Smoothing by convolution with g enables to consider again sets of particles rep- 
resented by delta-functions in p. 

To disprove the possibility of splitting W 2 into a sum of a positive and a 
positive definite function one may use the pu of equ. (IT51) embedded into M 2 , 

Pd(x,v) = p D (x)5(y). 

Due to the smoothing of W r by g and due to its decrease given by the decrease 
of h, the integral J W 2 Pd is finite: 

/oo +2 
W 2 p D (x) dx = W 2 (0) + 2W 2 {l)p{\) + W ^\ hu + n \)^ n ) 

v=\ n =-2 
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The bounded support of / and g is needed here as it was in equ. dHJ). The 
exponential decrease implies 



The "const." factor involves the integrals over f 2 and g 2 , the error term O(^) 
gives the difference between e~ 5eu /5u and e~ £ ( 5u+n '/(5v + n). The summations 
over v and n give 



The first part is negative and increases without limit when e — > 0, while the other 
term remains finite. So W% with small e can not be a sum of positive and positive 
definite functions. 

3 Mathematical background 

Only in applying Proposition[T]in statistical mechanics the Thermodynamic Limit 
is considered, not yet in the investigations of "stability". Moreover, in the re- 
formulation in [3] Definition there is no mentioning of "particles" . What is used 
of properties of space are: A distance relation between points and an invariant 
measure. This allows for a more general version of the definition, concerning 
functions on groups. We keep the notation we used above: x and y are elements 
of the group, their "group product" is x + y, the "inverse" of x is —x. 

5 DEFINITION. Consider a bounded continuous real valued function V on a 
locally compact abelian group G which has the Haar measure dx. V is stable, if 



for every finite positive Borel measure p(x)dx. 

Stable functions can be added, multiplied by positive numbers, and limits 
may be formed. So they form a closed convex cone, which we call STB. This 
cone STB contains POS, the cone of positive functions, also PDF, the cone of 
positive definite functions and sums thereof. 



An investigation of the relations between these cones may proceed via investi- 
gation of the dual cones (see [V641 IR621 IG03] ). The dual cones are subsets of 
V, the space of finite Borel measures n(x)dx, which is the dual space to V, the 






(15) 



STB d POS + PDF 



(16) 
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Banach space of bounded continuous functions. The dual cone to POS is POS', 
the set of finite positive Borel measures, dual to PDF is PDF', the set of finite 
positive definite Borel measures. The cone STB' is given as the closure of the 
cone of convex combinations of "correlation measures" 

H(x)= / p(x)p(y + x)dy, (17) 

J G 

i.e. convolutions of finite positive Borel measures p(x)dx with their reflected ver- 
sion p(—x)dx. These correlation measures are both positive and positive definite: 

STB' C POS' n PDF' (18) 

Now the question of equality or inequality in this relation is related to the 
central problem which is our concern in this investigation, the question of equality 
or inequality in (fl6j) . If the closed cone POS' fl PDF' contains an element pi which 
is not in the closed cone STB', then, by definition of "dual cone", there exists 
an element V G STB such that f V p, < 0, incompatible with a decompostion 
V=f + g, /ePOS, sgPDF. 

For the groups Z 2 , Z 3 , Z 4 there is equality in the equations (TIE]) and (TTB1) . but 
not for Z5. 

6 PROPOSITION. The intersection of POS' n PDF' with the plane 

2) . . .//(2))|/i(0) = 1} is completely characterized by its extremal points 
(0,0,1,0,0), (0,7,1,7,0), (7,0,1,0,7), (1,1,1,1,1), with 7 = (VE-l)/2 = 
l/(2|cos4vr/5|). 

Proof. By using Bochner's theorem and analyzing the Fourier transform 

2 

p,(k) = M^)e~ 27rfen/5 . (19) 

n=-2 

□ 

On the other hand there is a bound for STB' which cuts off a triangular subset 
of this convex quadrangle: 

7 LEMMA. Each element o/STB' obeys the inequality 

2 

Mi) < E ^V 4 - ( 2 °) 

n=-2 

Proof. STB' is defined by its extremal rays, formed as correlation measures of 
positive densities. 



p G STB', fj, extremal <^> 3p > 0, p,{n) = p(m)p(m + 



n 
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Assume, w.l.o.g., that p(— 1) > p(— 2). then 



p(l) = [p(-l) + p(l)] • [p(-2) + p(0) + p(2)] - [p(-l) - p(-2)] p(2) - p(-2)p(l) 



Remark: Also p(2) obeys this inequality and p(— 1) = p(l), p(— 2) = p(2). 
Closer inspection reveals moreover two rounded edges of STB'. 

Now the extremal point with p(n) as in equation (Q with v = is outside 
this boundary. And V{n) as in equation ([HD is an element of STB, but outside of 
POS+PDF. 

4 Conclusion 

For pair potentials which are bounded continuous functions the property of being 
"stable" can be reformulated without mention of particles. In this way it can be 
studied for abstract abelian groups. At the heart of the present investigation is 
the observation of a function V in Z5 which is stable, but indecomposable into a 
sum of positive and positive definite functions. This function V can also be used 
on Z. With some smoothing it can be used on R, and in damped periodically 
extended, rotationally symmetrized and again smoothed form on M 2 . Of course 
it is possible find sets of other examples nearby. So V(— 1) = V(l) in Theorem H] 
could be a little bit higher than —1. Only at — (y/E+ l)/4 « —0.8 does it become 
decomposable. 

The construction of a rotationally invariant example for dimension two is not 
so simple. A nicer one, or one for higher dimension, is not yet known. 
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Here s = P( m ) ' x 
Observe ^2 n p(n) = s 2 . 



, S \ / S \ s 

<(^-x)(- + x)<j. 
[p(-2)+p(0)+p(2)-p(-l)-p(l)]/2. 



□ 



